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Using kinetic theory for homogeneous collisionless magnetized plasmas, we present an extended 
review of the plasma waves and instabilities and discuss the anisotropic response of generalized 
relativistic dielectric tensor and Onsager symmetry properties for arbitrary distribution functions. 
In general, we observe that for such plasmas only those electromagnetic modes whose magnetic-field 
perturbations are perpendicular to the ambient magnetic field, i.e., Bi 1. Bo, are effected by 
the anisotropy. However, in oblique propagation all modes do show such anisotropic effects. 
Considering the non-relativistic bi-Maxwellian distribution and studying the relevant components 
of the general dielectric tensor under appropriate conditions, we derive the dispersion relations for 
various modes and instabilities. We show that only the electromagnetic R- and L- waves, those 
derived from them (i.e., the whistler mode, pure Alfven mode, firehose instability, and whistler 
instability), and the O-mode are affected by thermal anisotropies, since they satisfy the required 
condition Bi±Bo. By contrast, the perpendicularly propagating X-mode and the modes derived 
from it ( the pure transverse X-mode and Bernstein mode) show no such effect. In general, we note 
that the thermal anisotropy modifies the parallel propagating modes via the parallel acoustic effect, 
while it modifies the perpendicular propagating modes via the Larmor-radius effect. In oblique 
propagation for kinetic Alfven waves, the thermal anisotropy affects the kinetic regime more than 
it affects the inertial regime. The generalized fast mode exhibits two distinct acoustic effects, one 
in the direction parallel to the ambient magnetic field and the other in the direction perpendicular 
to it. In the fast-mode instability, the magneto-sonic wave causes suppression of the firehose 
instability. We discuss all these propagation characteristics and present graphic illustrations. The 
threshold conditions for different instabilities are also obtained. 
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I. INTRODUCTION 



Extensive studies have been conducted over the years using a wide variety of particle distributions to derive the 
general dielectric tensor so as to study various modes and instabilities [l|-[lH. Brambilla 0j also noted that, in a 
homogeneous collisionless plasma, the dielectric tensor always satisfies the Onsaeer symmetry relations independently 
of the particular form of the equilibrium distribution function. Gaelzer et al., [J| presented a detailed derivation of 
the effective longitudinal dielectric constant for plasmas in inhomogeneous magnetic fields. Ziebell and Schneider Q 
provided a general expression for an effective dielectric tensor satisfying the Onsager symmetry in an inhomogeneous 
plasma. Miyamoto |.6], choosing an anisotropic streaming Maxwellian plasma, obtained the general dielectric tensor 
and employed it to derive the dispersion relation for different instabilities. 

As is well known, the temperature is anisotropic in several environments, such as the solar wind, solar corona, auroral 
ionosphere, magnetosphere, astrophysical and space plasmas, and plasmas produced in recent experimental develop- 
ments. These naturally occurring and laboratory plasmas generally have bi-Maxwellian or nearly bi-Maxwellian 
velocity distributions. Noci et al., analyzed the data from the outer solar corona to show that the velocity dis- 
tribution function becomes anisotropic beyond 1.8 solar radii. Pagel et al., [l^ observed that the core population of 
solar wind is well described by a bi-Maxwellian distribution function, while the halo component is best modeled by a 
bi-Kappa distribution function. Schunk and Watkins 14] showed that the anisotropy in the electron-temperature dis- 
tribution in the ionosphere develops above 2500km and increases with altitude . Masood et al., [l^ analysed that the 
electron- velocity distribution in the Earth's magneto-sheath is thermally anisotropic. Cremaschini et al., [l6| showed 
that the accretion-disk plasmas around the compact and massive objects are characterized by thermally anisotropic 
velocity distribution functions . In the laboratory, the incident high-energy laser wave can make the resulting plasma 
thermally anisotropic, because the plasma is more intensely heated in the direction of the laser-wave electric field 

MM- 

The thermal anisotropy not only infiuences the propagation of various modes, but also induces various types of 
electromagnetic instabilities in a variety of plasmas including space and astrophysical plasma, fusion plasmas ( both 
magnetic and inertial confinement) as well as in the plasma created by highly intense free electron x-ray laser pulses. 
In particular, the stability analysis of whistler wave and Alfvenic modes (pure Alfven wave, magnetosonic wave, 
kinetic Alfven waves, inertial Alfven wave, etc.) have received special attention. Whistler waves are electromagnetic 
waves in magnetized plasmas with frequencies below the cyclotron frequency. Whistlers are naturally produced by 
lightning discharges in thunderstorms. When produced near the north or south pole, they can travel from one pole 
to another along the Earth's magnetic lines of force through the ionosphere. As a result, lightning hashes in the 
Southern Hemisphere can be observed in the Northern Hemisphere [l!|. Plasma- wave instruments on Voyager 2 have 
detected whistler-mode emissions inside and outside the magnetosphere of Saturn. Whistlers are also observed to 



propagate through self-created ducts in magnetospheres [22[. Whistler modes are used to induce radio-frequency 




plasma discharges, and to heat plasmas in tokamaks'23' and spheromaks[24|. The Weibel instability and the Weibel 
instability in an ambient magnetic field (i.e., the whistler instability), which have been known for several decades 
are of significant interest. Studies considering different particle-velocity distributions have relied on the Weibel 
(33} and whistler instabilities [33]-[i3 to explain physical processes in different plasma environments. 
Alfven waves are believed to play major roles in certain astrophysical processes in magnetized plasmas such as 
ones found in the environments of stars and interstellar clouds [ijj- Many phenomena in the solar atmosphere or 
heliosphere, planetary and cometary magnetospheres, cometary tails. Earth's ionosphere, etc., can also be regarded 
as manifestations of linear or nonlinear Alfven waves. Due to their incompressibility and low reflectivity in the solar 
atmosphere (44| . Alfven waves have been invoked as the most promising wave mechanism to explain the heating of Sun's 
outer atmosphere, or corona, to millions of degrees and the acceleration of the solar wind to hundreds of kilometers 
per second[4^. Gekelman has studied Alfven waves and their relationship to space observations in the laboratory 
[46| . One of the most important instabilities that excite Alfven waves in a hot plasma is the so-called fire-hose 
instability, driven by the plasma temperature anisotropy. Since the thermal anisotropy is an intrinsic characteristic 
of magnetized plasmas, especially in the cases of collisionless plasmas such as astrophysical and space plasmas, the 
fire-hose instability has many implications in these plasmas, and has been discussed in the literature on the basis of 
the Vlasov kinetic theory. The fire-hose instability in magnetized thermal plasmas has been studied by a number 
of authors 47|-[50|. Kinetic Alfven waves (KAWs) are produced due to charge separation when the perpendicular 



wavelength becomes comparable to the ion gyroradius and show dispersive character for oblique propagation j51| 
Bashir et al., discussed the effect of thermal anisotropy on the propagation characteristics of KAWs in both the 



kinetic and the inerial regimes[52l|. Various other aspects of the KAWs have been studied by several authors [53|-[70|. 
Recent theoretical and experimental advances focused on the wave propagation of KAWs in different regimes have 
associated the morphology of the wave with the generation mechanism proposed by Gekelman et al., [7l|. 

Here, we present an extended review on plasma waves and instabilities and discuss the anisotropic response of 
the general relativistic dielectric tensor as well as the Onsager symmetric properties for a homogeneous magnetized 
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collisionless plasma for an arbitrary distribution function. Considering non-relativistic bi-Maxwellian distributions, 
we simplify the analytical expressions for the components of dielectric tensor in the limit |^„q| 3> 1, where £^na — 
(a; — n^a) / {k\\Vt\\a)] this allows easier switching to parallel or perpendicular propagation. We moreover obtain the 
dispersion relations for various modes and instabilities, along with their graphical representation, from the simplified 
components of the dielectric tensor under the appropriate conditions. For example, (i) for parallel propagation, 
we derive the dispersion relations for the R- and L-waves, whistler wave, Alfven wave, Langmuir wave, Alfven wave 
instability, whistler instabilities, and Weibel instabilities, along with the conditions for instability; (ii) for perpendicular 
propagation, the general dispersion relations for the X-mode, 0-mode and Bernstein mode are derived; (iii) for oblique 
propagation, we derive the general dispersion relations for the Kinetic Alfven waves (KAWs) in the kinetic and inertial 
regimes, and for the fast mode. In addition, we recover a number of special cases for the general fast-mode dispersion 
relation under appropriate conditions. We also discuss the fast mode instability for oblique propagation. 

The plan of the paper is as follows. In Section II, we present the general relativistic dielectric tensor and then 
employ the non-relativistic bi-Maxwellian distribution function to simplify the general dispersion relation for various 
modes and instabilities under appropriate conditions, along with their graphical representations. In Section III we 
briefly summarize and discuss the results. 



II. MATHEMATICAL MODEL 



We start out with the relativistic Vlasov equation 

dfa , dfa 



VXB^ dfa, 



dt Ox ' ' c 

where the relativistic momentum p and velocity v are related by 

p = -ymv 7 = I 1 ^ 
cp 

V = = 

Equation (1), combined with the Maxwell Equations 

VxE = — — 
c ot 

„ „ IdE 47r, 
VxB = -— + — J, 

c ot c 

describe the dynamics of the plasma system. 



dp 



(1) 



(2) 

(3) 
(4) 



Linearizing Eqs. ([T}, ([3]), and taking the Fourier-Laplace transform, assuming the ambient magnetic field Bp 
along z- direction and the wavenumber vector k in the x-z plane, we obtain the perturbed distribution function 



and the dyadic equation 



{(s + (0 — 0') + ik±v±{s'm(j> — sine/)')} 



(s^ + c^k^) Ei-c2k(k.Ei) + 47r.sJ = 



dp' 



(5) 



Given the distribution function fia, we can calculate the current density 



where ^ is the conductivity tensor. 
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We can rewrite the dyadic Eq. ([5]) in terms of the dielectric permittivity tensor ^[6| as 



Ej = 0, 



where 



On H = Oi. 



(6) 



(7) 



with 



Stt f f 



Xi 



^^;_L \Jn{z)\ 



in 



z 



w±J„(z)j;(z) 



J„(z)J^(z) 



V 



the only nonzero component of Lij behig 



U|| / dfoa dfoa 



v± 



dpi 



v± 



d 



P\\ 



Here we have used the foUowing properties of the Bessel functions: 

oo oo oo 

n— — oc n— — oo 

and defined 

xi = i ^r: — ^ ^\\- ''^-L 



n— — oo 



dpj_ 



dp\\ 



N 



ck 



is the total refractive index, whose components are Nij — chi^ j/^, while uj'i^ = 47rnoQ,e2/ rriQa is the 
non-relativistic plasma frequency, and 57^ [— rt^aM — QchBq / {•'ymQac)] is the relativistic cyclotron frequency. 

Equation ([7]) gives the general relativistic dielectric tensor in a homogeneous magnetized plasma for an arbitrary 
equilibrium distribution function /qq. From Eq.(l7]), it is clear that xi ( only for fcii 7^ 0) and L^z contain the anisotropic 



modification factor 



^11 



dfoa 



dfoa 



This factor stems from the term (v x Bi) 



dfaa 



in the linearized Vlasov 



dpj_ dpii ) ap 
equation, which vanishes for isotropic velocity distributions or under electrostatic conditions, i.e., for k||E. The 
same term however survives for those electromagnetic modes in which the ambient magnetic field is perpendicular to 
the perturbed magnetic field, i.e., for Bi _L Bp, for any arbitrary anisotropic distribution (from non-relativistic to 
relativistic regimes). From the Maxwell Equations, it is clear that generation of perturbed magnetic field depends 
upon the propagation direction and the polarizations of electric field. As illustrated by Fig. 1 the R- and L-waves 
and the waves derived from them for parallel propagation, and the 0-mode for perpendicular propagation satisfy 
the requirement that the ambient magnetic field be perpendicular to the perturbed magnetic field. The same is 
evident from Eq.Q, since xi exhibits the thermal anisotropy effect for ^ only. In the case of the 0-mode (i.e., 
Ezz — -/V^ = 0), the anisotropy effect comes from the L^z component. For parallel propagation of electrostatic modes 
(i.e., for tzz = 0), a contribution from the anisotropic effect is found in and another is found in L^, but the two 
contributions cancel each other. The obliquely propagating electromagnetic modes always contain the anisotropic 
modification. 

In the following section we will prove that for a non-relativistic bi-Maxwellian distribution, only the modes whose 
perturbed magnetic field is perpendicular to the ambient field are modified by the thermal anisotropy. 
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A remarkable property of the general dielectric tensor is described by the Onsager symmetry relations, 

*^xy ^yx^ ^xz ^zx ^^nd 6^^ 

which is always satisfied by the hot-plasma dielectric tensor, independently of the choice of equilibrium distribution 
function /qq,. 

The above discussion on the anisotropic response of dielectric tensor and Onsager symmetry relations are valid for 
arbitrary anisotropic distributions and cover magnetized-plasma environments ranging from the non-relativistic to 
the ultra-relativistic regimes. 



A. General Dielectric Tensor for Non-relativistic Bi-Maxwellian Distribution 



We shall now derive the general dielectric tensor for a non-relativistic bi-Maxwellian plasma. The non-relativistic 
bi-Maxwellian distribution function is 



foa — 



1 



— exp 



Pi 



P 



(8) 



Once the integrations over p|| and p± are carried out, with the help of the Onsager symmetry relations, the following 
general dispersion relation is obtained 52 1 



^xx 



xy 



^xy ^xz 

^yy ~ ^yz 



TV? 



= 



(9) 



where 



-^-Y."-stf: Sr..(A,){^..zK,„)-q^ 1)} (10) 



n — — Qo 



Q n — — oo ^ ^ ^ V II' 



~-yy 



i-E#^o" E r„(A„)e„„z'(e„„)|i + ^(^^-i)^ (12) 

71 — — 00 



2 



— — oo 



sxz^-Y.^^^-^ E ^^|eo. + e.J^~i)^^'(^„.) (14) 

^ ijj^ Vt^^ ^ V2Aa t \ ^l|a / J 

^^^ = *E?^^ E \/5r;(A.)|eo. + 6.o(|^-i)}^'(Cn.) (15) 

Q,*-La ^j^_„^ L \||q /J 

Here, eij is the permittivity dielectric tensor. The parameters in Eqs. (fTQ)) -(fT5l) are defined as follows: 



Cna = , A„ = '^*ll" = (2T||Jm„)'/' and v,^^ = (2r^„/m„) (16) 

'<^||^t||a ^^'Oq 

The results of the integration overpx are expressed in terms of the functions r„(AQ,) and V^(\a)^ which are related to 
the modified Bessel function (A „) by the expressions r„(AQ) = e""^" ln{)^a) andF^Aa) = e~'*'° (/'^(A^) — /^(Aq)). 
The integration over p|| introduces the plasma dispersion functions Z {£,na) and Z' {£,na) [z3- For the isotropic case, 
the dielectric permittivity tensor eij reduces to the usual textbook form [l|, Q . 
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B. Dispersion Relations for various Modes and Instabilities 



In the limit \£,na \ S> 1, the general dispersion relation is further simplified, and the dispersion relations for various 
modes and instabilities can be derived. In that limit, the plasma dispersion function assumes the form 



Z {S,na) — ~7 



5 

net 



E 



(2^ + 1)!! / fc||Wt 



{21 + 1) 2' \uj - jin, 



21 + 1 



for|e„a| > 1 



For ICnctl ^ 1, the components of the tensor in Eq.® take the form 



2 CO 9 



{21 + 1) (T^. 



{2l + l)\\ { w(/c||Vt||J' 



2/+1 



^ (2^ + 1) \{u:-nn^^) 



£yy — 1 



^2 ^ 



2 oo 



2i+i 
n'r„(Aa) 



•||"*||c 



{uj + niloa) 

; s 2/+2N 



, (2^ + 1) 



2i+i 

,2 



1^*11 



21 



'(^P*llJ' 



\2/+l 



\ 2/+2 / X 2/+2N 



{UJ + niloa 



,2; + ! 



a; — nSloa , 



UJ - 



T, 



Til 



... = 1 - E %ro(A„)eo^„ z' (^0.) - E — E rn(Ao) E 



,2 oo 



J=0 



(2^ + 1)!! 
2' 



(fc||Wt„ 



v2; 



21+1 



{uj + nQoa 



,21 + 1 



(^iPtiiJ 



21 



(fc||Wt,| 



{uj - nfioa) 



2/+1 



21+1 



UJ 



exy- Syx- ,„2 Z^^^nl^aJZ^ CO/ 4- ^ ^ . ^2/+l 



v2! 



,21 



+ 



71 = 1 

{21 + 1) (T+^ 



2'+i V T\ 



{21 + 1) {{uj-nn^X^^ {uj + nVL^^) 



21 + 1 



1 



Oct 



w - 



Oa 



(17) 



(18) 



(19) 



(20) 



E' 



^ nr„(A„)^ (2/ + 1)!! 



^ ""iic ,fr^i V2Aa 



i=0 



21 



(fc||Vt|| 



, 21+1 



2;+i 



(oj — nQoa^ 



21+2 



{uj + nf^Oa) 



(w - nf^oa) 



2Z+2 



,21 + 1 



21+2 



(w + nf^oQ 



2/+1 



,2i+2 



(21) 



7 



2 oo 



Vt± LU 



n = l 



2 

21+1 



2' 

X 2i + l 



Tllry 



, 2J+1 



.21 + 1 



{uj - nfloa) 



21+2 



2;+2 



(22) 

In Eqs. (IT7|) - (I^^ . we have separated the n = terms from those with n > and n < 0. For n = 0, we have retained 
the plasma dispersion function to be used later. We have also used the symmetry properties of the modified Bessel 
function, i .e., r„(AQ) = r_„(AQ.) and r^(AQ.) = r'_„(AQ,). Each component of the tensor contributes to oblique 
propagation, since both components of the wave-number vector, and fc|| are present via the r„(Act) and Z{£^na) 
functions, respectively, and contain the effects of thermal anisotropy. 

Sections III B 11 III B 21 and IIB 3 discuss the parallel, perpendicular and oblique propagations, respectively, for the 
limiting case |^„q,| ^ 1. 



1. Parallel Propagation 
Taking fcx = in Eqs. ([T7 | - ([22)) . we obtain the equations 



1 



/ .2 °° 



(2Z + 1)!! I ^{h^^y , ^(fc|Pt„J 



.21 



[21 + 1 ) (T+ 
2i+i 



I" 



LU — rioo 



2;+2N 



-xy 



__.^<.^(2/ + l)!! / w(fc||«t||J 



2/+1 



(2^ + 1) /7±a _ ^^ / 



2'+i \ Til 



2i+2N 



and Exz = Sy 







Here we have expanded the modified Bessel function as follows 



1 



An 



and rQ(AQ) 



-1 



and noted that only the n = I terms survive for k+ = 0. 

Therefore, for parallel propagation, the general dispersion relation in Eq.® reduces to the form 



e.z = 0, 



which decouples the modes 



from the electromagnetic modes 



0, 



-^11 — £xx i i^xy 



(23) 

(24) 
(25) 



(26) 

(27) 
(28) 
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The upper (+) and the lower (— ) signs correspond to the R- and L- waves, respectively. While the dispersion rela- 
tion for R- and L- waves retain the thermal anisotropy effect, the electrostatic mode exhibits no such effect, as expected. 

a. R- and L- Waves 

Using Eqs. ([23l) and ([24|) . we may rewrite the dispersion relation for the R- and L- waves in Eq. ((28l) in the form 

>p<f. (2/ + l)!! / fcp«,i [ u: (21 + 1) fT^^ \ ( \] 

^ (2^ + 1) Uif^oc.; \^±r!oa 2'+i VT||„ J \{uj±noo.rj f ^ ' 

This equation represents the general dispersion relation for R- and L- waves incorporating the thermal anisotropy 
as well as the higher-order thermal effects parallel to the ambient magnetic field. To retrieve the cold plasma result, 
we may assume the plasma to be isotropic and take only the leading term, i.e., I — 0. Further, with i7oQ = , the 
structure of both the R- and the L- waves vanishes altogether and we obtain the dispersion relation for an unmagnetized 
plasma. Since the dielectric tensor for more than one species is additive, we can also solve the dispersion relation for 
a multi-species plasma. 

The graphical representations of electron R- and L- waves with / = in Eq. (30) are shown in Figs. 2 and 3, 
respectively. The dotted curve represents the standard isotropic R- and L- wave and the other curves depict deviation 
from the isotropic case. We observe that the phase velocity of the R- and L- waves increases as the magnitude 
of anisotropy is increased, but the cutoff point remains unaffected. For the R-wave, the resonance disappears as 
anisotropy is switched on. 

h. Whistler Mode 

If we keep only the leading term (i.e., Z = 0) in Eq. ip^ and assume that cj^ ^ c^fc|, the dispersion relation for 
an electron R-wave reduces to 



2 V Tiu UJ 



(30) 



where 



A| = ^ = ^^,with T/j^-^. (31) 

In the low-frequency range f2oi < uj < floe , Eq. pop describes the general whistler mode 



^pe I V^lle J 2 




Til ' / 



(32) 



UJ 



Here we have retained the first-order term in . 

For longer wavelengths, i.e., uj^^ 3> c^fc^, this dispersion relation reduces to the form 



^pe r ' V^lle 'J 2 



(33) 



This is the dispersion relation for a whistler mode including the electronic thermal anisotropy effect. If we neglect 

UJ 

the second order term in the denominator on the right-hand side, which is due to the first order term in — — , we 

recover the results of Lazar et al., [42|. Figure 4 shows how the phase speed of the whistler wave is enhanced as the 
thermal anisotropy increases. 
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c. Non-Resonant Whistler Instability (\(,na \ ^ IJ 

By considering only the leading term (i.e., / = 0) in the R-wave and assuming that ^ c^fc^ in Eq. ([29l) . we obtain 
the following expressions for the real and imaginary parts of lo: 



2 (c2fc2+a.2^ 



(34) 



and 



2(^-1) fc|i;2^ (c2fc2+^2^)-f72^a.2^ 



2 (c2A:|+a.2^ 



(35) 



in agreement with Lee et al., [37| and Lazar et al., [42|. 

The instability occurs for wave numbers satisfying the condition 



,2,2 _ ^ 



02 



(36) 



From Fig. 5, we see that the growth rate increases as the thermal anisotropy increases, while the instability threshold 
is lowered. For any specific anisotropy value, the growth rate initially rises for small /c|| and then saturates at large 

/C||. 

From Eq. I|35p . we see that the ambient magnetic field suppresses the growth rate. In the limit of large thermal 
anisotropy, i.e., for ^ 3> 1, the dispersion relation reduces to 



Oe'^pe 



2 (c2fc2+,,2^ 



(37) 



If the magnetic field is neglected, Eq. (|37| reduces to the simpler expression 



V2 



c2fc^+ILj2^ 



which is the well-known result due to Weibel [25| . 
d. Alfven Wave 



We now consider the dispersion relation for a two-component plasma consisting of electrons and singly charged 
ions. In the low- frequency regime lj < Qi, we keep only the leading term, i.e., Z = 0, to rewrite Eq. (|29p in the form 



'4 



7±c 

Til 



1 ± 



2uj 



Tn 



2uj 



(38) 



Here we have assumed that the two species are at the same parallel temperature and taken advantage of the small 
mass ratio, rrie/mi <Si 1. 

In the isotropic case, Eq. (p8| reduces to the standard pure Alfven mode. For frequencies satisfying ut ^ iloi, the 
quadratic structure in oj vanishes, and both the R- and the L- waves follow the same dispersion relation. 



uj = fell Vk 



Til 



- 1 



T 
'Tw 



- 1 



(39) 
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a result derived by Schlickeiser and Skoda [49j . 

In Eqs. ([38]) and (pQ]) . the Alfven wave frequency is modified by acoustic effects due to the thermal anisotropy of the 
electrons and ions. The Alfven phase velocity can be enhanced or reduced, depending on the intensity and signatures 
of the anisotropy. 

From Eq. ((39|) . it is evident that the constraint 



2 ( T_ie 



1 - 



T 



Til 



makes the Alfven mode unstable, for then 



W = i fc|| Vyl 



1 - 



7±e 



1 - 



Til 



- 1 



(40) 



which describes the so-called fire-hose instability. The Alfven- wave frequency is real for T_Leji > T\\e,ii under which 
condition fire-hose instabilities cannot arise. 



e. Electrostatic Waves 
From Eq. (j27p . i.e., ezz = 0, we obtain the dispersion relation for Langmuir waves 

2 2 I 7 2 2 



(41) 



where we have used the condition |^oe| ^ 1 ioi expanding the plasma dispersion function Z{^o 
For Vt\\i <^uj /k\\ <^ Ut||e and kjXj^^ ^ 1, the equation e^z — yields the ion-acoustic mode 



2 I 2 2 



where we have used that 



1 T 

Z' (Coe) = -2 for ICoel « 1, Z' iU) = TT fo'^ l^o^l » 1 and cl = 



(42) 



The electrostatic modes, which are only affected by the temperature in the direction of propagation, are hence 
insensitive to the thermal anisotropy. 



2. Perpendicular Propagation 

Letting fc|| — and noting that only the leading terms in the Z-summation (i.e., I — 0) survive in Eqs. (|17p - ([22)) . we 
obtain the expressions 



n = 1 



^ zz 



1 - E ^ E - r;.,A.) + 2 E r;,A,) ,44, 

1 - E ^r„(A„) - 2$: £ r„(A.) (l + ^ j (45) 



a n—1 



Q n — 1 





In the expressions for eyy and ezz, Eqs. (IT51) and respectively, we have expanded the plasma-dispersion functions 
Z(eoa) and Z'iCoc) for l^oal > 1- 
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Thus the general dispersion relation in Eq. ([9]) reduces to 



{e-zz 



Ni 







(47) 



The above equation shows two decoupled modes: the Extraordinary (X-) mode 



(48) 



and the Ordinary (0-) mode 



Ni = tzz 



(49) 



If we assume the non-diagonal components txy to be much smaller than the diagonal components t^x and ^yy, the 
X-mode dispersion relation (P5|) yields two decoupled modes: the purely transversal X-mode 

(50) 



Ni 



and the Bernstein mode 



0. 



(51) 



Here, like in parallel propagation, we observe that the thermal anisotropy affects the dispersion relation of the 
0-mode [Ea. (|i^ ]. but does not affect the X-mode [Ea. (|I5)) ] or the Bernstein mode [Eg.dST])]. 
We will next discuss these three modes in more detail. 

a. X- Mode 



With the help of Eqs. (HU and we rewrite the X-mode dispersion relation, Eq. (pS]) in the form 



,21,2 



c'k 



n = 1 



n"r„(Ac,) 



2Aq, V'^JyXa^ 



2o2 

Oa 



(^2 „2 

E fr„(A„) 



7)202 



(52) 



which is the general form of the X-mode. 

Figure 6 depicts the behavior of X-mode for an electron plasma. Panels (a), (c), and (e) show the behaviors for 
n — 1, and the harmonics n = 2, and 3, respectively, with Ae = 0.05, while panels (b), (d), and (f) correspond to 
rt = 1, 2, and 3, respectively, with Ae = 8.0. The number of cutoffs and resonances increases with n. The distinction 
between the resonance and the cutoff points is more prominent in the right panels, i.e., for the large argument. The 
cutoff and resonance points shift towards lower frequencies and hence expands the propagation domain as Ae grows 
from 0.05 to 8.0. 



The pure transverse X-mode (i.e.,c^fc^/a;^ = ^yy)i given by Eq. (|50p . may be written as 



= 1 - 



7^ 2-^ 
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(53) 



The two X-mode dispersion relations, Eqs. (15^ and ([55]) . contain the standard modified Bessel function. We 
can solve these relations numerically to describe the X-mode in detail. Alternatively, we can seek an approximate 
analytical solution results by expanding the Bessel function for small argument, which converts Eqs. (j52[) . and (j53p 
to the form 
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(54) 
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and 



, ,2 <x> 2/1 \ 2n — 2 / 2 \ 

2^ ,„2 2^ niy 2^0. J W - n^^L J ' ^ ^ 



1 



respectively. 

Equation (|55p is the expression derived by Zaheer and Murtaza for pure transverse X-modes [TSj. Equation ([5l)) 
with the sums over n on the right-hand side truncated at n ~ 1 reduces to the textbook expression for the dispersion 
relation for the general X-mode in electron plasmas (see, e. g., Chen [l]) : 



(56) 



h. Bernstein Wave 

In view of the expression for Eq. (|33]), the dispersion relation for the Bernstein wave t^x — in Eq. ([51]) takes 
the form 



a ^ ' n — 1 ^ ^ 

Figure 7 depicts solutions of this equation. The dotted line represents the solution for Ag = 1 and the other three 
curves represent Ag = 0.5, 1.5, and 2.0. The resonance points are fixed, but the cutoff points shift from right to 
left and reduce the propagation domain as Ae grows. At higher harmonics (not shown), the cutoff points become 
independent of Ae. 

c. O- Mode 

Given the expression for e^z, Eq. (|45|) . the dispersion relation for the 0-mode in Ea. (l49l) is given by the equality 

a Q n — 1 



Again we expand the Bessel function for small argument and rewrite this equation in the form 



^2 ^Z^ ^2 i 2rin« y lw2_„2f^2 Y Ti „ / W2 - ^2^(2 f ^^^i 



n! V 2fioa / l^^^-"^f^0Q V Tj-Q / w2 - n2ri2^ 



This is a general dispersion relation for the 0-mode with higher order thermal effects, including the consequences of 
thermal anisotropy. 



3. Oblique Propagation 



In the low-frequency, long parallel- wavelength regime, the non-diagonal components of the tensor Sij are negligibly 
small. Therefore the dispersion relation in Eq. (9), for the kinetic Alfven waves (KAWs) can be written as 



0. 



In this regime, the fast mode is described by the expression 



'VV 



- = 0, 



(60) 



(61) 



and is hence decoupled from the KAWs. 

The electric field vector E and the wavenumber vector k are coplanar for the KAWs [68l|-[70l|. 

Reference [s^l presents a detailed derivation of the expressions describing the oblique propagation of the KAWs in 
the kinetic and inertial limits. Here we present a simpler derivation. To highlight the effects of thermal anisotropy 
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and to show the wave frequency w as a function of the perpendicular and parallel wave vectors and we present 
a 3-D graphical representation of the simplified results. To the best of our knowledge, our results for the general fast 
mode and the fast mode instability are new. 

General expressions for the components of the tensor Sij were presented in Section III Al We substitute Eq. pUj) for 
Exx and Eq. (fT2)) for Ezz and consider the low- frequency limit w ^ fi^ , k^^Xjj^ ^ 1, ^ to reduce Eq. (l60l) to 
the equality 



ro(Ae) ^, , ro(A,) 




(62) 



with the shorthand 



Here, the anisotropy terms ipi and -02 are 

^2 
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2 ^ t\\a 
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and 



Til, \ , _ /T, 



^-2 = 1^ (1 - ro(A,)) - 1 j + (1 - ro(A,)) - ij I (64) 

Although in Eq. (j62p , we have included contributions from both species for terms containing anisotropic effect but 
we have neglected the electronic contribution from the isotropic part of the e^x component by using the small mass 
ratio, rrie/mi <C 1. 

Next, we shall derive the dispersion relation for the kinetic Alfven wave in the kinetic (i.e., Vt\\i uj /k\\ <^ Vt\\e 
and {rrie/mi) <^ P — {c^s\J^a) ^ 1) ^"^^ the incrtial (wt||e,i <SC /^|| , P ^ nie/rrii) limits. The kinetic (inertial) 
limit is defined with reference to the dominant thermal (inertial) electronic effect. 

a. KAWs in the Kinetic Limit 

To discuss the kinetic limit of KAWs, we assume the parallel phase velocity of the wave to be less than the parallel 
thermal velocity of the electrons, yet greater than the parallel thermal velocity of the ions, i.e., <^ to /k^^ <C v^^^ , 
and the plasma to have low /3, i.e., {nie/mi) <C /3 = {cl^^ )/Va) ^ 1- Under these conditions, the dispersion relation in 
Eq. ([62)1 reduces to 



^7.2 2\M-Ie ,\ , T||j / 3 2 2 



where 
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In deriving Ea. (|55)) we have assumed the gyro-radii to be small, so that 

3 A 1 3 

To{Xe,^) Kl- Xe.i + -Xii and ——2^—— = - ^ ^ c± 1 -|- -fclp^.i- (67) 

Equation ([62]) describes the general kinetic Alfven wave, in which ip'i measures the deviation from the isotropic 
limit. Illustrative plots are presented in Fig. 8. 
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b. KAWs in the Inertial Limit 



To obtain the inertial limit of KAWs, we assume that the parallel phase velocity of the wave is greater than the parallel 
thermal velocities of both the electrons and the ions, i.e., Wt||e,j ^ /^|| , P ^ {me/'^i), and that the gyroradii are 
small. Under these constraints, the dispersion relation in Eq. ((62)) yields the following expression for the modified 
kinetic Alfven wave in the inertial regime: 



The inertial Alfven wave is modified by the acoustic effect, in turn caused by the thermal anisotropy. Since 
j3 <C me/mi only extraordinarily large thermal anisotropics will show appreciable effects. 



Til, 



T\\e 



(68) 



c. Fast Mode 



We take advantage of the low- frequency limit u <^ J7oa and of the Bessel function identity 2 ^ r„(AQ) = 1 — ro(AQ) 

71=1 

and keep only the leading term (/ = 0) in the sum over I on the right-hand side of Eq. (j20p . The fast mode dispersion 
relation in Eq. (|6T|) ( i.e., Syy — N"^ = 0) then reduces to 



n= 1 



2 




+2^-i^A„r(,(A„)|l+(^-l) ^. (69) 

Here we have expanded the plasma dispersion function Z (^oa) for |^oa| S> 1. 

To further simplify Eq. , we let <C and take the small gyro radius limit (in which the argument of the 
Bessel function is small), so that the term containing r^(AQ) becomes negligible. This leads to the result 

where we have used that rQ(AQ) = —1 and defined = 'i;^j^^/(2r2oQ). 

For an electron-ion plasma, this dispersion relation can be written in the form 

= k^Vl + kl^ + k^cl + (71) 



where 

Equation ([nT) . the modified dispersion relation for the fast mode, incorporates the thermal anisotropy effects of the 
two species. The parallel- propagating term displays an acoustic effect associated with the parallel temperature, i.e., 
which is sensitive to the thermal anisotropy. By contrast, the acoustic effect c^^ associated with the perpendicular- 
propagating term depends on the perpendicular temperature, but receives no contribution from the temperature 
anisotropy. In oblique propagation, the general fast mode has therefore two distinct acoustic effects, one in the 
parallel and the other in the perpendicular direction. The enhancement or the reduction of the fast mode frequency 
depends upon the strength and signature of the thermal anisotropics of the two species. 

Figure 9 represents the general fast mode for small and for large /3|| 's. The anisotropy effect becomes more prominent 
when is large. 

Three instances of the general fast-mode dispersion relation ([71]) deserve special mention: 

(i) For parallel propagation (i.e., k± =0), Eq. ([7T]) gives the Alfven- wave dispersion relation in Eq. (|39p . namely, 



(72) 
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Here we have assumed that the two species have the same parallel temperature. 

(ii) For perpendicular propagation (fc|| — 0), Eq. ([7T|) yields the standard dispersion relation for the magneto-sonic 
mode, 

(iii) In the isotropic limit i.e., T±a = T\\a = Eq. (|71l) reduces to the standard fast mode dispersion relation 

= evi + kici 

where 

Cs = 



d. Fast Mode Instability 
Under the constraint 
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(73) 



the fast mode, Eq. (f7T|) . becomes unstable, since 

r,2 



Ti 



1 - 



T 

Til 



- 1 



VI 



(74) 



The inequality requires that /3|| = c^^/Vj > 1 and that T±a < T||q,. 

In Eq. (1741) , the parallel-propagating part represents the Alfven wave with thermal anisotropy while the perpendic- 
ular part represents the magneto-sonic wave, which tends to suppress instabilities. Therefore, fire-hose instabilities 
may be damped for small perpendicular wavelengths, as demonstrated by Fig. 10. 



In the other limit, |^„q| < 1, we expand the plasma dispersion function as 



4£4 
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(2; + l)!! V 
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and proceed to discuss the parallel-propagating whistler instability. 

e. Resonant Whistler Instability ( \S,na \ l£ ^) 
With Eqs. ((TU)) . pT|) . and (IT5|) . the dispersion relation for R- waves in electron plasma becomes 



2;+i 



(75) 
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(76) 



where we have kept only the leading term of the plasma dispersion function for 



< 1. 

Under the subluminal condition, i.e., for a; <C cfc , the real and imaginary parts of w are 
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(77) 



and 
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This defines the whistler instability, which occurs for waves satisfying the wavenumber condition 



,2 . ^pe fT^e 



1 



) 



The magnetic field generates real oscillations, but the growth rate remains unaffected. These results coincide with 
those due to Lazar et al., [12] ■ Figure 11 exhibits the effect of anisotropy on the whistler instability. For the field-free 
case, i.e., Bq = 0, only the purely growing Weibel instability arises j25|]. 



On the basis of kinetic theory, we have extensively reviewed plasma waves and instabilities. In particular, we find 
that for any anisotropic equilibrium distribution, here included the non-relativistic, relativistic, and ultra-relativistic 
magnetized collisionless homogeneous plasmas, the electrostatic modes are insensitive to thermal anisotropics, which 
affect only the electromagnetic modes with magnetic field perturbations perpendicular to the ambient magnetic field. 

We have also derived expressions describing a general dielectric tensor for a magnetized non-relativistic bi- 
Maxwellian plasma. In alignment with the salient physical features of the problem, our analysis showed that 
the dielectric tensor receives separate contributions from the integrations over the perpendicular and the parallel 
momentum components. The integration over the perpendicular components leads to the functions r„(AQ) and 
T'^^Xa), which are related to the modified Bessel function /„(Aq), while the integration over the parallel components 
introduces the plasma dispersion functions Z (^na) and Z' (^„a). From the resulting general dielectric tensor we have 
obtained the dispersion relations for a variety of modes and instabilities. 

We have also expanded the plasma dispersion functions in the limit l^nal 3> 1 to derive relatively simple analytical 
expressions for the components of dielectric tensor and to find the dispersion relations for those modes and instabilities. 
For example, (i) for parallel propagation, we have expanded the modified Bessel function as fcj^ — >■ and note that 
only n = 1 terms survive. We have also derived the dispersion relations for the R- and L- waves, whistler wave, 
Alfven wave, Langmuir wave, non-resonant whistler instability, Weibel instability and Alfven-wave instability; (ii) for 
perpendicular propagation, we have let fc|| — > and noted that only the leading term {I = 0) in the sum over / survives 
to derive the general dispersion relations for the X-, 0-, and Bernstein modes, (iii) For oblique propagation, we have 
found the general dispersion relations for the kinetic Alfven wave in the kinetic and the inertial regimes and for the 
fast mode. Additionally, we have discussed a few special cases of the fast mode and the fast mode instability. For 
\Cna\ < 1, we have derived the parallel- propagating whistler instability. 

In both limiting cases, as expected, the electrostatic modes are unaffected by the thermal anisotropy. Among 
the electromagnetic modes, the parallel-propagating modes, such as the R- and L- waves and the modes derived 
from them (the whistler mode, pure Alfven mode, firehose instability, and Weibel instability) are affected, while the 
perpendicularly-propagating modes, such as the X-mode and the modes derived from it (the pure transverse X- and 
Bernstein modes) are not. The 0-mode, however, is sensitive to thermal anisotropics. 

The thermal anisotropy affects the parallel propagating modes via the acoustic effect, while it affects the perpen- 
dicular propagating modes via Larmor-radius effects. For oblique propagation, both effects appear, additively, in the 
Alfvenic modes. 

The effects of thermal anisotropics are more prominent in the kinetic limit than the inertial limit. The parallel- 
propagating term of the fast mode displays an acoustic effect with parallel temperature, i.e., controlled by the 
thermal anisotropy. On the other hand, the perpendicular-propagating term displays an acoustic effect with perpen- 
dicular temperature, i.e., c^^, but it is independent of the thermal anisotropy. In oblique propagation, the general 
fast mode presents two distinct acoustic effects, in the parallel and perpendicular directions. 

For both the resonant and the non-resonant cases, the whistler and Weibel instabilities can only exist in environments 
with T±a > T\\a whereas the firehose and the general fast-mode instabilities can only arise for Tuq, > T±a- 

In conclusion, we find that the thermal anisotropy affects only the modes with magnetic perturbations perpendicular 
to the ambient magnetic field i.e., with Bi _L Bq. The anisotropy can either enhance or reduce the frequency domain 
of waves and instabilities, depending on its strength and signature. 

Our results may prove useful for studies of the thermally anisotropic environments frequently found in astrophysical, 
space, and even laboratory plasmas. 
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Captions: 

Figure 1. This figure shows the geometrical representation of different modes. In all panels, the z axis is aligned 
with the ambient field Bq- The arrows indicate the directions of the electric field E\ , magnetic field B\ (blue bold 
arrow) and the wave vector fc (red bold arrow). Panels (a) and (b) show an electromagnetic mode and an electrostatic 
mode propagating in the z direction and. Panel (c) depicts an extraordinary mode, the coupling of an electrostatic 
and an electromagnetic modes, propagating in the x direction. Panel (d) shows an ordinary mode propagating in the 
X direction. 

Figure 2. This is the graphical illustration of circularly polarized electron R-wave. The normalized phase veloc- 
ity w^/c^/c^ as function of normalized wave frequency w/fioe which depicts the deviation from isotropic case with 
the change of thermal anisotropy values ( i.e., A=0(dottcd), 5(grccn), 20(blue) , 40(red)) choosing the parameters 
Wj,e/^^Oe = l-2,i't||/c = 0.2 & A = (Tj_e/T||e) — 1. It shows that the phase velocity of R-wave increases with the 
increase in temperature anisotropy. 

Figure 3. This is the graphical illustration of circularly polarized electron L-wavc. The normalized phase veloc- 
ity w^/c^A:^ as function of normalized wave frequency w/f2oe which depicts the deviation from isotropic case with 
the change of thermal anisotropy values ( i.e., A=0(dotted), 5(green), 20(blue) , 40(red)) choosing the parameters 
Wpe/f^oe = 1.2,wt||/c = 0.2 & j4 = (Tj_e/T||e) — 1. It shows that the phase velocity of L-wave increases with the 
increase in temperature anisotropy. 

Figure 4. This figure illustrate the graph of the normalized wave frequency w/^e of the whistler mode as a function 
of normalized wavcmimber cfc/wpe with the parameter ojpe/f^Oe = lOjWfu/c = 0.05 & A = (Tie/T||e) — 1 by choosing 
different values of thermal anisotropy ( i.e., A=0 (dotted), 5 (green), 20 (blue) , 40 (red)). It shows that the wave 
frequency increases with the increase in temperature anisotropy. 

Figure 5. The normalized growth rate Im w/cjpe of non-resonant whistler instability (|^„q| ^ 1) as a fimction of 
normalized wavcnumbcr cfc/wpe with the parameter Wpe/^^oe = l-0,i't||/c — 0.1 A = (Tie/T'iie) ~ 1 by choosing 
different values of thermal anisotropy ( i.e., A= (Black), 5 (green), 10 (red), 20 (Blue)). By increasing the temper- 
ature anisotropy. The growth rate increases and threshold on wave number shifts towards the longer wavelength and 
thus enlarges the wave vector domain. This instability exist only in the plasma environments where with Tlq > T\\a- 

Figure 6. This figure depicts the behavior of X-mode for an electron plasma. Panels (a), (c), and (e) show the 
behaviors for n = 1 (dotted), and the harmonics n = 2 (green), andn = 3 (red) , respectively, with Ae = 0.05, while 
panels (b), (d), and (f) correspond to n = 1, 2, and 3, respectively, with Ae = 8.0. The number of cutoffs and 
resonances increases with n. The distinction between the resonance and the cutoff points is more prominent in the 
right panels, i.e., for the large argument. The cutoff and resonance points shift towards lower frequencies and hence 
expand the propagation domain as Ag grows from 0.05 to 8.0. 

Figure 7. The electron Bernstein mode is illustrated graphically. The dotted line represents the solution for Ae = 1 
and the other three curves represent Ae = 0.5 (black), 1.0 (green), 1.5 (blue), and 2.0 (red). The resonance points 
are fixed, but the cutoff points shift from right to left and reduce the propagation domain as Ae grows. At higher 
harmonics (not shown), the cutoff points become independent of A,.. 

Figure 8. The normalized frequency of KAWs in kinetic limit u) (fc||,fcj_) /wpi as a function of perpendicular and 
parallel normalized wavenumber i.e., k_\_pi/ujpi), ~ 0.05 & Ag = Ai = A — (Tj_/r||) — 1 by choosing different 
values of thermal anisotropy ( i.e., A=0 (red), 10 (blue), 20(green) , 40(purple)). The frequency increases with the 
increase of temperature anisotropy value and the effect of temperature anisotropy is more prominent in high beta 
plasma environments than the low beta one where T±a > 2^||c«- 

Figure 9. The normalized fast mode frequency co [k\\ , k±) /cOpi as a funtion of perpendicular and parallel normalized 
wavenumber i.e., (fc||Uyi/ ojpi & k±VA/^pi) with the parameters A^ = Ai = A = (Ti/Th) — 1; (a)/3|| = 0.05 & {b)(3\\ 
= 10.0 by choosing different values of thermal anisotropy ( i.e., A=0 (red), 10 (blue), 15(green) , 30 purple). The 
frequency increases with the increase of temperature anisotropy value and the effect of temperature anisotropy is more 
prominent in high beta plasma environments than the low beta one where T±a > 

Figure 10. The normalized growth rate Imw (fc||,fc_L) /ujpi of fast mode instability as a function of perpendicular 
and parallel normalized wavenumber i.e., (fc||Uyi/ & k±VA/i^pi) with the parameters A^. = Ai = A ^ (T^^/T±) & 
/3|| = 3.0 by choosing different values of thermal anisotropy ( i.e., A=l (red), 2 (blue), 4 (green) , 9 (purple)). The 
growth rate enhances with the increase of temperature anisotropy but the perpendicular wave vector stabilizes this 
instability i.e., the magneto-sonic wave suppresses the fire-hose instability which exist only in the high beta plasma 
environment where with T±a > T\\a- 

Figure 11. The normalized growth rate Im w/wpe of non-resonant whistler instability {\^na\ < 1) as a function of 
normalized wavenumber ck/upe with the parameters Vt^\/c = 0.1 & A = (T±e/T\\e) — 1 by choosing different values of 
thermal anisotropy i.e., A= (Black), 5 (green), 10 (red), 20 (Blue). By increasing the temperature anisotropy. The 
growth rate increases and threshold on wave number shifts towards the longer wavelength and thus enlarges the wave 
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vector domain. This instability exist only in the plasma environment where with T±a > ^||a 
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